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1. Introduction 

The study of quasi-one dimensional strongly correlated electron systems (SCES's) 
has received a lot of attention in the past few years. This interest was mainly trig- 
gered by the synthesis of materials which in a wide range of temperatures can be 
well modeled by a 3D system of (almost) decoupled spin chains, spin ladders or more 
generally, Hubbard chains and ladders [1]. An important reason that has put these 
ID systems again on the scene is the appearance of the so-called "stripe" phases 
which have shown up in very different contexts, such as high-T c cuprates [2], integer 
quantum Hall effect systems at high Landau levels [3], materials showing colossal 
magnetoresistance (manganites) [4], etc. 

The present article contains a brief review of theoretical work on magnetization 
plateaux in SCES's in ID. The presentation reflects the authors' perspective and due 
to lack of space, not all important contributions to this subject can be mentioned. The 
first studies on this issue were performed by Hida [5] and Okamoto [6] in an attempt 
to describe some organic compounds with periodic couplings. After a few other 
isolated cases were studied [7,8], Oshikawa and collaborators [9] have undertaken 
the first systematic study of this problem and, by extending the Lieb-Schultz-Mattis 
theorem to systems in a magnetic field, they provided a necessary condition for the 
appearance of magnetization plateaux in ID systems. When the magnetization (M) 
is normalized to saturation values ±1, this condition for the appearance of a plateau 
with magnetization (M) can be cast in the form 

SV (1 — (M)) G Z . (1) 
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Here S is the size of the local spin and V the number of spins in the unit cell for 
the translation operator acting on the magnetization (M) groundstate 1 . It should be 
noted that translational invariance can be spontaneously broken in the groundstate 
and then V would be larger than the unit cell of the Hamiltonian. 

Spin ladders in a magnetic field constitute a class of systems where the full 
phase diagram was explored and where it was checked when the necessary condition 
(1) becomes also sufficient [15-20]. 

A famous result on pure ID spin chains is Haldane's conjecture [21] (see also 
[22] for a new field theoretical approach). It states that such chains with integer 
S are gapful while those with half-integer S remain gapless, corresponding to the 
condition (1) with (M) = and V = 1 since a spin gap is equivalent to an (M) = 
plateau. This feature is indeed confirmed in seminal numerical studies of Heisenberg 
chains [23,24], but no other plateaux are found for these systems. Other non-trivial 
plateaux, which for S > 1 would be permitted by (1), arise only when the model is 
modified, e.g. by adding a single-ion anisotropy [25,26]. 

The existence of real materials with trimer constituents [27-31] motivated the 
investigation of the magnetization process of trimerized [5,6,32], frustrated trimer- 
ized [33,34], quadrumerized [35] and more generally periodically modulated S = 1/2 
Heisenberg spin chains with period p (so-called p-merized chains) [36-40]. Ladders 
consisting of both staggered and non-staggered dimerization along the chains were 
also studied [41]. 

From the experimental point of view, one of the most exciting materials is 
NH4CUCI3 where just two plateaux with (M) = 1/4 and 3/4 have been observed [42]. 
The crystal structure of NH4CUCI3 suggests to model it as a two-leg zig-zag ladder, 
however with room for some modifications. While one can indeed obtain plateaux 
with (M) = 1/4 and 3/4 in a zig-zag ladder with dimerized legs [43,44], neither this 
point of view nor other proposals [45] have so far resulted in a really satisfactory 
theoretical description of the experimental observations. In short, if NH4CUCI3 is 
a quasi-one-dimensional system, eq. (1) should apply with 5 = 1/2 since the spin 
is carried by Cu 2+ ions. Then one needs V = 8 in order to permit plateaux with 
(M) = 1/4 and 3/4. However, V = 8 would permit also plateaux with (M) = and 
1/2 which are not observed experimentally [42], but are very pronounced e.g. in the 
dimerized zig-zag ladder [44]. 

When charge degrees of freedom are dynamical, as in Hubbard chains, the sit- 
uation turns out to be quite different and very interesting, although the effect of a 
magnetic field has been studied so far essentially only in two cases: An integrable 
spin-S 1 generalization of the t — J model [46] and doped Hubbard chains with peri- 
odically modulated hopping matrix elements or on-site energies [47]. In addition to 
fully gapped situations (both a charge gap and a spin gap), situations arise in which 
magnetization plateaux appear at irrational doping-dependent values of the mag- 
netization, in contrast to what we have just reviewed for purely magnetic systems. 
Interestingly, under such circumstances, gapless degrees of freedom are present and 
consequently certain superconducting correlations have power law behavior [47]. 



1 In cases where the sites carry different spins [10-14], the combination SV should 
be replaced by the maximal spin in the unit cell. 
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Recently, the effect of disorder on plateau systems, such as the p-merized XX Z 
chains, has also been studied. It was found that, for binary disorder, plateaux do 
not disappear but are rather shifted in a precise amount which can be predicted by 
means of a simple argument to be [48] (M) = 1 + | (p — 1), where p is the strength 
of disorder. Conversely, continuous distributions of disorder erase completely the 
plateau structure. Interestingly, the effect of disorder on the susceptibility leads to 
an even-odd effect similar to what is found in iV-leg ladders [49] . 



2. Field theory approach 

The field theory description of spin chains is a useful technique for treating 
problems of weakly coupled spin ladders, weak dimerization in chains, etc. as we will 
illustrate below for the case of S = 1/2. For higher spin chains, non-Abelian bosoniza- 
tion has proven to be better suited. For spin 1/2 chains, the Abelian bosonization 
technique has proven to be very efficient. It describes the low energy, large scale 
behavior of the system, and can be extended to the case of an easy axis anisotropy. 
More specifically, the continuum limit of the Hamiltonian 



Hxxz = J X <AS X S X+1 + - (S+S x+1 + S x S+ +1 ) | - h £| (2) 

x=l { ' x=l 

is given by the Tomonaga-Luttinger Hamiltonian 

H=±Jdx (vK(dJ) 2 + ^{d^f) . (3) 

The bosonic field (j) 1 and its dual 4> l are given by the sum and difference of the lightcone 
components, respectively. The constant K = K((M),A) governs the conformal 
dimensions of the bosonic vertex operators and can be obtained exactly from the 
Bethe Ansatz solution of the XX Z chain (see e.g. [16] for a detailed summary). One 
has K = 1 for the SU(2) symmetric case (A = 1) and it is related to the radius R of 
[16] by K~ x = 2nR 2 . 

In terms of these fields, the spin operators read 

S^ = ^=d x (p + a:cos(2k F x + V2^(P):+^p-, (4) 
V 2n 2 



: (-l) x : e ±iv ^ (bcos(2k F x + v 7 ^) + c) : , (5) 

where the colons denote normal ordering with respect to the groundstate with magne- 
tization (M). The Fermi momentum kp is related to the magnetization of the chain 
as kp = (1 — (M))tt/2. An XX Z anisotropy and/or the external magnetic field mod- 
ify the scaling dimensions of the physical fields through K and the commensurability 
properties of the spin operators, as can be seen from (4), (5). 

If one introduces periodically modulated couplings between the spins J x = J if 
x 7^ np and J x = (1 — S) J = J' if x = np, this amounts to the following perturbation 
for the original Hamiltonian in the continuum limit [36] 

H per t = X J dx cos(2pkpx + V2~TV(/)) , (6) 
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where A is proportional to 5 J. 

Another interesting situation is the one of spin ladders [15,16], where iV identical 
chains are coupled with a transversal coupling J' 

N a=N 
a=l x,a=l 

For simplicity we used here periodic boundary conditions (PBC's) along the trans- 
verse direction. One obtains in the continuum a collection of identical Hamiltonians 
like (3), with perturbation terms which couple the fields of the different chains. After 
a careful renormalization group (RG) analysis, one can show that at most one degree 
of freedom, given by the combination of fields (pp = S a 0a, remains massless. The 
large scale effective action for the ladder systems is then given again by a Hamiltonian 
(3) for (j)o and the perturbation term 

H per t = \ J dx cos(2Nk F x + V2ncf) D ) , (7) 

where kp = (1 — (M))ir/2 is related to the total magnetization (M). 

In both cases, the key point if to identify the values of the magnetization for 
which the perturbation operators (6) or (7) can play an important role. In fact, this 
operator is commensurate at values of the magnetization given by (1) with S = 1/2 
and V = p or V = N, respectively. If this operator turns out to be also relevant in 
the RG sense (this depends on the parameters of the effective Hamiltonian (3), the 
model will have a finite gap, implying a plateau in the magnetization curve. The 
generalization to p-merization in A^-leg ladders results in V = pN, though in some 
cases the topology of couplings has to be carefully analyzed [41] 

The case of a single p-merized chain can be generalized to the presence of charge 
carriers [47]. Let us consider the Hubbard model, representing interacting electrons 
with spin 1/2. The Hamiltonian is given by 



x 

h 



x.a x 
x,a x 



where t(x) and (i(x) are taken as periodic in the variable x with period p. The 
bosonized version of this Hamiltonian can be written as a Gaussian part, given by 



J2jjdx [(d^) 2 + (dj t ) 2 ] 



l = C,S 



where the fields labeled by c and s are particular combinations of the bosonic fields 
dictated by the Bethe Ansatz solution of the homogeneous model. The most relevant 
of several perturbation terms is given by 

Opert = Ai sin[— +pnnx — \/tt^4>c] cosfv^yr^s] + A2 sin[7m + 2pnnx — V4n^(p c ] . (9) 
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We have written the perturbation term for the case of zero magnetic field for simplic- 
ity. Using the same arguments of commensurability and relevance as above, we can 
show that for pn G Z, we have a charge gap. If the condition is further constrained 
to pn/2 G Z, for zero magnetic field we have also a spin gap, implying a plateau with 
(M) = 0. In the case of non-zero magnetic field, if one of the conditions 

|(n±(M))GZ, 

is satisfied, and the doping is kept fixed (as is natural from the point of view of ex- 
perimental realizations of doped systems), the system has a magnetization plateau, 
but still exhibits massless behavior as well, e.g. in the specific heat which vanishes 
linearly as the temperature goes to zero. If both conditions are simultaneously sat- 
isfied, the system is gapped in the charge and spin sectors; this situation is in fact 
the generalization to arbitrary doping of the results for p-merized Heisenberg chains 
discussed above. 



3. Strong coupling arguments 

The magnetization process is easy to understand if the system decouples into 
clusters of V sites [15,16,37,47]. These 'strongly coupled' clusters magnetize indepen- 
dently such that at zero temperature the magnetization (M) can only take finitely 
many values. For a spin-system with spin S they are subject to the quantization con- 
dition (1) where V is directly the number of spins in one of the decoupled clusters. 
Using (higher order) series expansions around this decoupling point [15,16,34,37], the 
quantization condition obtained in this simple manner can be argued to be valid also 
for more general parameters. 

In its simplest form, this argument explains magnetization plateaux which arise 
because the unit cell of the Hamiltonian contains several spins. However, the ar- 
gument can be refined to account also for spontaneous breaking of translational 
symmetry by a period of two. The transitions between plateau states can be treated 
by degenerate perturbation theory, leading to effective Hamiltonians which in many 
cases turn out to be effectively an XX Z chain (2) if only the lowest order(s) are kept 
[50-54,17,43,34]. If the effective XX Z anisotropy turns out to be sufficiently large 
(A > 1), a gap opens and translational symmetry is spontaneously broken. In this 
manner one finds a further plateau precisely in the middle between the two values of 
(M) predicted by considering just the decoupling limit itself. This illustrates that V 
in (1) should be taken from the unit cell of the groundstate whose size is an integer 
multiple (in general larger than 1) of the unit cell of the Hamiltonian. 

Amusingly, first-order strong-coupling results for the magnetization process re- 
main exact over a finite range of parameters in certain systems with local conservation 
laws [12,55]. 

4. Numerical approaches 

In many cases coupling constants do not lie in any of the weak- or strong- 
coupling regimes discussed before. In such cases, accurate results can be obtained 
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Figure 1. Magnetization curves of (a) the planar 3-leg ladder at J' = 3 J 
[15] and (b) the quadrumerized chain with J' = J/2 [36]. The diamonds 
show results of (a) fourth- [16] and (b) second-order [37] strong-coupling 
expansions for the plateau boundaries. The bold line is an extrapolation 
to the thermodynamic limit. 



only numerically. Since on the one hand, vector spaces are large already for moderate 
system sizes and on the other hand, only a few extreme eigenvalues are needed to 
compute the zero-temperature magnetization curves, one frequently resorts to the 
Lanczos method as was already done in pioneering works [23,24]. In the present 
situation the magnetic field couples to a conserved quantity S* ot . One therefore only 
needs the groundstate energy E(S* ot , h = 0) at h = in each of the magnetization 
subspaces with S^ ot G {0, 1, ... , SL} (here L is the total number of spins). Then one 
can readily obtain the energy in a finite field h through the relation E(S* ot ,h) = 
E(S* ot ,Q) — hS^ ot and then construct the groundstate magnetization curve. Note 
that the magnetization has steps on a finite system since the possible values of the 
magnetization are quantized. 

Among the systems whose magnetization process has been studied by the Lanc- 
zos method, the two-leg zig-zag ladder (which is equivalent to the J\-Ji chain) has 
a particularly long tradition [56-61]. However, plateaux with (M) ^ are ob- 
served only for more than two legs [43,38] or if the coupling constants are modulated 
[41,44,50,61,62]. 

Fig. la) shows the magnetization curve of an S = 1/2 3-leg ladder with open 
boundary conditions along the rungs [15]. As expected, there is a clear plateau 
with (M) = 1/3 and the fourth-order strong coupling series for the boundaries of 
the plateau [16] are in good agreement with the Lanczos results. The effect of 
p-merization is illustrated in Fig. lb) which shows the magnetization curve of a 
quadrumerized S = 1/2 Heisenberg chain [36] . Here one observes two clear plateaux 
at (M) = and 1/2, as expected. The series for their boundaries [37] are in good 
agreement with the numerical data even if they are only of second order. Moreover, 
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Figure 2. Magnetization curves of (a) the S = 1 Heisenberg chain and 
(b) the 5 = 3/2-1/2 ferrimagnetic Heisenberg chain. The steps show the 
curves obtained for a finite system with L = 60 sites by DMRG [12] while 
the smooth curves denote an extrapolation to the thermodynamic limit. 
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in the latter case bosonization predicts the plateaux to open with power laws in S 
which can indeed be verified numerically [36]. 

An altenative method is DMRG [63] which allows to study larger systems. The 
number of applications of this method or variants thereof to the magnetization process 
is increasing steadily [11,12,17-19,33,38,62,64-67]. The only technical complication is 
that the infinite system algorithm does not give good approximations to incommen- 
surate structures and therefore the finite system algorithm should be used if one 
wants to compute the complete magnetization curve. Fig. 2 illustrates this method 
with two examples [12]: (a) The magnetization curve of the S = 1 Heisenberg chain 
exhibits only an (M) = plateau, corresponding to the Haldane gap [21], and (b) 
the magnetization curve of the S = 3/2-1/2 ferrimagnetic Heisenberg chain which 
exhibits a spontaneous magnetization and thus a plateau with (M) = 1/2. 

Long XX chains in a transversal field can be treated by the Jordan- Wigner 
transformation (see e.g. [40]). This is particularly useful for systems with quenched 
disorder [48]. Fig. 3 illustrates the effect of a disordered distribution of couplings 
with a p-periodically modulated background on the magnetic behavior. Here, we 
show the magnetization curves of a quadrumerized XX chain under different disorder 
strengths p of a binary distribution P(Ji) = pS(Ji — J') + (1 — p)S(Ji — J — 7^J), 
where 7i = 7 , (—7) ii i = pn , (i ^ pn). One observes the appearance of a rather 
nontrivial phenomenon, namely, the shift of the magnetization values (not necessarily 
rational), for which certain plateaux emerge, as compared to the pure system. One 
can show by means of a simple real space decimation procedure [48] that these new 
plateaux appear at (M) = 1 + 2(p — l)/p for XX Z chains which in the special case 
of XX chains is consistent with the numerical results. Continuous distributions for 
the disorder wipe out completely the plateau structure. 
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Figure 3. Magnetization curves of quadrumerized XX chains with 5 x 10 4 
sites in disordered binary backgrounds of strength p = 0.2, 0.4, 0.6, 0.8 (solid 
lines in ascending order) [48]. The left and rightmost dotted lines denote 
respectively the pure uniform and pure modulated cases p = 1 and p = 0. 



5. Experimental realizations 

Two concrete plateau substances have already been mentioned earlier: Clear 
plateaux with (M) = 1/4 and 3/4 have been observed in NH4CUCI3 [42], but in 
this case a thorough theoretical understanding is still lacking. Conversely, a plateau 
with (M) = 1/3 is theoretically expected [33,34] in Cu 3 Cl 6 (H20)2-2H 8 C 4 S02, but 
magnetization experiments [31] would have to be extended to slightly higher fields 
to confirm its presence. One example where experimental observations do actually 
confirm theoretical predictions [7,8] are the plateaux with (M) = 1/2 and which 
have been observed in the bond-alternating S = 1 chain compound [Ni 2 (Medpt)2(/U- 
ox)(//-N 3 )]C1O 4 -0.5H 2 O [65]. 

iVdeg spin ladders are realized in modifications of high-T c materials [1] . However, 
magnetization plateaux are predicted for these systems at fields of (several) thousand 
Tesla which is clearly outside the presently accessible experimental range. One al- 
ternative is provided by organic compounds and indeed magnetization experiments 
have been successfully performed on the two-leg ladder material Cu2(C5Hi 2 N 2 )2Cl4 
(see e.g. [52]). Vanadium oxides may in general provide another alternative and e.g. 
some PBC 3- leg ladder is suspected in the recently synthesized compound Na2Vs07 
[68]. However, there is already one experimental observation which can be interpreted 
in terms of a 3- leg ladder: A plateau-like feature at (M) = 1/3 has been observed 
in CsCuCl3 for a magnetic field perpendicular to the crystal axis [69]. Although 
CsCuCl3 is strictly speaking three-dimensional, this plateau can be understood in the 
approximation of a weakly anisotropic PBC 3-leg ladder with ferromagnetic rungs 
[67]. 

In general, it is intriguing that doping can push magnetization plateaux into 
the low-field region at least for p-merized Hubbard chains [47] . This points towards 
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the possibility of observing magnetization plateaux in doped systems whose parent 
compounds would exhibit plateaux only at inaccessibly large fields. 
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